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INTRODUCTION 


Nearly  two  decades  of  magnetospheric  research  have  established 
that  charged  particles  from  the  solar  wind  are  energized  inside  the  magneto- 
sphere, although  the  mechanism  of  their  energization  is  not  yet  precisely 
determined.  Hot  magnetospheric  plasma  are  often  observed  (DeForest  and 
Mcllwain,  1971;  and  references  therein)  in  the  equatorial  region  on  auroral 
field  lines.  Observations  of  energetic  electrons  at  low  altitudes  on  the  same 
auroral  field  lines  (Evans,  1975;  auid  references  therein)  and,  more  recently, 
observations  of  energetic  ions  of  ionospheric  orgin  (Hultqvist  et  al,  , 1971; 
Shelley  et  al.  , 1976;  Mizera  et  al.  , 1976)  confirm  that  magnetospheric  and 
ionospheric  plasmas  interact  strongly  in  the  auroral  region.  In  particular, 
beam-like  pitch-angle  and  energy  distributions  are  often  observed  and  seem 
to  indicate  that  such  energetic  particles  have  been  accelerated  or  retarded 
by  an  electric  field  E parallel  to  the  magnetic  field  B . Moreover,  it  has 
been  pointed  out  (Evans,  1975)  that  the  characteristics  of  such  energetic 
electron  events  may  perhaps  be  roughly  classified  into  two  categories: 

(a)  transient  events  for  which  the  flux  distributions  indicate  diffusive  energy 
gain,  and  (b)  quiescent  quasi- static  events  for  which  the  pitch- angle  and  flux 
distributions  show  a beam-like  behavior.  Since  the  magnetospheric  plasma 
during  a substorm  event  is  highly  disturbed,  the  existence  of  the  first  cate- 
gory of  event  is  not  at  all  surprising;  the  diffusive  energy  gain  might  be 
associated  with  the  anomalous  resistance  that  can  occur  along  a magnetic 
field  line  as  consequence  of  several  instability  mechanisms  (e.g.,  Perkins, 
1968;  Kindel  and  Kennel,  1971  ; Papadopoulos  and  Coffey,  1974).  The 
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second  category  of  event,  which  does  not  show  diffusive  energy  gain,  seems 
to  require  a quasi- static  electric  potential  difference  of  up  to  several  kilo- 
volts between  the  ionosphere  and  the  equator  in  order  to  account  for  the  beam- 
like behavior  of  electron  flux  distributions.  Although  a parallel  electric 
field  can  arise  as  the  result  of  mapping  perpendicular  electric  fields  down- 
ward into  the  ionosphere,  where  collisional  resistance  along  and  collisional 
conductance  across  the  magnetic  field  may  become  important  (Chiu,  1974), 
such  a process  is  limited  to  the  lower  ionosphere  and  cannot  account  for 
the  beam-like  characteristics  of  energetic- electron  precipitation.  In  regions 
of  open  field  lines,  charge- separation  effects  related  to  the  polar  wind  also 
give  rise  to  a parallel  electric  field  (Banks  and  Holzer,  1968);  however,  the 
magnitude  of  such  a total  potential  drop  is  far  less  than  the  potential  differ- 
ences of  up  to  several  kilovolts  that  are  inferred  to  occur  between  the  iono- 
sphere and  the  equator  along  certain  (presumably  closed)  auroral  field  lines. 
Thus,  the  evident  existence  of  quasi- static  parallel  electric  fields  in  the 
auroral  region  remains  a puzzle  in  auroral  physics.  In  this  paper  we  pro- 
pose to  consider  the  origin  of  such  parallel  electrostatic  fields  and  their 
mutual  consistency  with  auroral  plasma. 

Basically,  there  are  two  mechanisms  by  which  a quasi-static  electric 
potential  difference  V can  be  established  along  a magnetic  field  line  in  a 
collisionless  plasma.  On  the  one  hand,  one  can  appeal  to  quasi- stable  double 
layers  that  have  been  produced  in  bounded  laboratory  plasmas  (Quon  and  Wong, 
1976)  in  which  there  are  interpenetrating  plasma  streams  involving  ions  and 
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electrons  reflected  from  the  walls.  Although  it  has  been  conjectured  (Block, 

1975)  that  such  double  layers  may  also  exist  in  the  magnetosphere,  where  plasma 
boundaries  are  somewhat  amorphous,  we  shall  temporarily  defer  considera- 
tion of  the  double-layer  mechanism  until  its  existence  and  stability  can  be 
verified  for  magnetospheric  plasma.  On  the  other  hand,  Alfven  and  Faltham- 
mar  (1963,  pp.  163-167)  have  pointed  out  that  an  anisotropic  collisionless 
plasma  in  a magnetic  field  can  be  in  quasi- neutral  equilibrium  without  a 
parallel  electric  field  only  if  the  magnetic  field  is  homogeneous  or  if  the 
pitch- angle  anisotropy  is  the  same  for  both  electrons  and  ions.  Otherwise, 
in  the  case  of  a dipolar  magnetic  field  such  as  the  earth's,  the  absence  of 
a parallel  electric  field  v/ould  result  in  different  distributions  of  ions  and 
electrons  along  the  field  line  if  the  equatorial  anisotropies  of  the  particles 
were  different.  Conversely,  charge  neutrality  of  the  plasma  along  the 
magnetic  field  line  demands  that  an  electrostatic  field  be  established  parallel 
to  the  magnetic  field  line.  The  electrostatic  potential  energy  difference 
|e|V^  between  the  ionosphere  (where  B = ) and  the  equator  (where  B = Bq  ) 

is  estimated  to  be  of  the  order  of  the  mirror  ratio  B^/B^  multiplied  by  the 
I mean  perpendicular  particle  energy  (Alfven  and  Falthammar,  1963;  Persson, 

j 1966)  , i.  e. , to  be  of  the  order  of  hundreds  of  keV  for  plasma  in  an  auroral 

flux  tube.  However,  the  magnetospheric- ionospheric  flux  tube  along  an 
auroral  field  line  contains  not  only  anisotropic  magnetospheric  plasma,  but 
also  thermal  plasma  extracted  or  evaporated  from  the  ionosphere  and  back- 

; scattered  electrons  of  intermediate  energy.  Since  the  various  components 

[ 
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of  colder  plasma  (here  assumed  isotropic  at  the  foot  of  the  field  line)  are 
expected  to  participate  in  the  maintenauice  of  charge  neutrality,  we  would 
expect  the  above  estimates  of  electrostatic  potential  difference  along  an 
auroral  field  line  (Alfven  and  Falthammar,  1963;  Persson,  1966)  to  be 
unrealistically  large.  Indeed,  calculations  of  the  parallel  electric  field  by 
application  of  the  principle  of  quasi- neutrality  (Lemaire  and  Scherer,  1973) 
yielded  electrostatic  potential  differences  of  only  a few  volts  between  the 
ionosphere  and  the  equator,  although  the  effects  of  backscattered  electrons 
and  of  pitch- angle  anisotropy  for  magnetospheric  plasma  were  ignored  in 
those  calculations.  In  view  of  the  importance  of  the  question,  we  have 
undertaken  to  re-examine  the  conditions  which  define  quasi- neutrality  for 
auroral  plasma,  especially  for  flux  tubes  in  the  high  latitude  trough  region  of  the 
ionosphere,  poleward  of  the  plasmapause,  where  the  energetic  electron 
observations  have  been  made. 

In  this  paper,  we  apply  the  principle  of  quasi-neutrality  to  calculate 
the  mutually  consistent  electrostatic  potential  and  particle  distributions 
along  an  auroral  field  line  which  is  populated  by  collisionless  anisotropic 
magnetospheric  plasma  and  by  plasma  extracted  or  backscattered  from  the 
ionosphere.  It  is  shown  'hat  a potential  difference  of  up  to  several  kilovolts 
between  the  equator  and  the  -ionosphere  may  thus  be  maintained  along  an 
auroral  field  line.  While  the  parallel  electric  field  has  the  proper  signature 
to  account  for  electron  precipitation  characteristics  (Evans,  1975),  it  also  may 
self- consistently  account  for  the  presence  of  ions  with  keV  energies 

at  the  upper  reaches  of  auroral  field  lines  (Shelley  et  al.  . 1976). 


Even  though  our  calculation  in  eimilar  to  that  of  Lemaire  and  Scherer 
(1973,  1977),  there  are  a number  of  i,rucial  differences  which  combine  to 
yield  a larger  potential  drop  along  the  field  line  in  our  work  than  in  theirs. 
Among  these,  three  essential  factors  peculiar  to  auroral  field  lines  in  the 
trough  region  are  most  important  in  accounting  for  the  differences;  indeed 
these  may  also  explain  why  the  beam- like  characteristics  of  energetic 
particles  are  not  observed  elsewhere.  First,  as  has  been  pointed  out  by 
Alfven  and  Faltheunmar  (1963),  different  anisotropies  for  electrons  and  pro- 
tons are  crucial  to  the  maintenance  of  a large  potential  drop  along  the  field 
line.  Since  ring-current  particles  injected  onto  the  auroral  field  line  are 
expected  to  be  anisotropic  in  pitch  angle,  we  preserve  such  a feature  in  our 
calculation.  Second,  copious  backscattered  electrons  are  observed  and  must 
be  considered  in  a quasi- neutrality  calculation.  Third,  the  number  density 
of  thermal  ions  at  the  ionosphere  is  an  important  boundary  condition  on  the 
problem,  since  the  ions  are  accelerated  upward  by  the  same  parallel  elec- 
tric field  that  accounts  for  downward  acceleration  of  electrons.  Since  the 
trough  region  of  the  ionosphere  is  considerably  depleted  relative  to  the 
average  ionosphere,  it  is  expected  that  the  effect  of  the  ionospheric  plasma, 
which  has  a tendency  to  reduce  the  potential  difference  along  a magnetic 
field  line,  is  correspondingly  minimized. 

We  shall  show  that  in  our  model  the  energetic  electron  flux  and 
pitch-angle  distributions  at  the  foot  of  the  field  line  are  essentially  similar 
to  those  found  in  the  previous  model  (Evans,  1975)  for  electron  precipi- 
tation undet  the  influence  of  arbitrarily  postulated  parallel  electric  fields. 


In  our  calculation,  however,  we  require  that  the  electrostatic  potential 
be  a solution  to  the  quasi-neUtrality  equation,  and  this  requirement  is 
found  to  be  an  important  constraint  on  the  admissibility  of  the  electro- 
static potential,  otherwise  postulated,  for  the  problem.  It  should  be 
noted  that  our  approach,  though  yielding  similar  flux  distribuf . ^ ns 
at  the  foot  of  the  field  line,  is  basically  different  from  Evans'  a'pro.ich. 
In  Evans’  model,  the  beam-like  characteristics  of  precipitating  electrons 
are  attribute<i  to  the  collimation  of  an  isotropic  electron  population 
by  an  arbitrarily  imposed  parallel  electric  field.  In  our  model,  the 
parallel  electric  field  arises  naturally  as  a consequence  of  the  quasi- 
neutrality requirement  imposed  on  the  various  particle  populations,  some 
of  which  are  assumed  to  be  anisotropic  at  the  equator.  Thus,  the 
mapping  relationships  between  electron  anisotropy  and  parallel  electric 
field  for  the  two  models  are  different. 
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QUASI-  NEU  TRALIT  Y 


The  principle  of  quae i- neutrality  is  commonly  employed  in  problems 

of  the  present  type  (e.g.  , Lemaire  and  Scherer,  1973).  One  considers 

a dipolar  magnetic  flux  tube  in  which  the  magnitude  of  B varies  monoton- 

ically  from  Bq  at  the  equator  to  B^  at  the  foot,  which  is  located  at  an 

appropriate  altitude  to  be  discussed  below.  We  construct  a kinetic  model 

in  which  the  ion  and  electron  distribution  functions,  f ^ and  f , are 

expressed  analytically  in  terms  of  the  velocity  components  Vj^  and  v^ 

relative  to  B , in  terms  of  the  local  electrostatic  potential  V s V(s), 

and  in  terms  of  the  local  magnetic  intensity  B s B(s),  where  s is  the 

8 

coordinate  that  measures  arc  length  of  the  field  line  from  the  equator.  The 
principle  of  quasi- neutrality  asserts  that 

[ t.(v,  . i v^)  . t_(v,,  V^)  ] dv,  = 0. 

(1) 

For  our  model,  in  which  and  f are  explicitly  constructed,  we  inter- 
pret (1)  as  an  equation  for  the  value  of  at  any  poi^t  s along  the  field 

line  ( 0 s 8 £ / ) under  the  convention  that  V vanishes  at  the  equator  (i.  e. , 
Vq  = 0 ).  Persson  (1966)  has  shown  that  a necessary  and  sufficient  con- 
dition for  V to  vanish  at  all  points  along  the  field  line  is  for  f and  f 
to  have  the  same  anisotropy,  i.  e.  , the  same  equatorial  pitch- angle  dis- 
tribution. Since  ring-current  electrons  and  protons  are  believed  to  be 
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scattered  in  pitch  angle  by  different  wave  modes  (Kennel  and  Petschok, 

1966;  Cornwall,  1966),  we  expect  that  their  pitch-angle  distributions  are 

not  the  same.  Thus,  an  electrostatic  potential  V ^ 0 is  required  in 

6 

order  to  satisfy  the  quasi-neutrality  condition  expressed  by  (1). 

Since  our  model  is  intended  for  application  in  the  topside  iono- 
sphere and  the  magnetosphere,  the  plasma  is  assumed  to  be  collision- 
less in  the  interval  0<  s s X.  This  assumption  of  collisionless  plasma 
in  the  flux  tube  entails  at  least  two  consequences  which  require  con- 
sideration at  the  outset.  On  the  one  hand,  the  distribution  functions 
f^  and  f are  required  to  be  mapped  along  the  field  line  in  accordance 
with  Liouville's  theorem,  the  detailed  consequences  of  which  will  be 
discussed  below.  On  the  other  hand,  we  are  obligated  to  select  X so 
that  all  the  charged-particle  populations  are  indeed  collisionless  throughout 
the  interval  0 5 s s X.  For  the  magnetospher ic  proton  population,  which 
is  obviously  collisionless  in  the  equatorial  region  of  the  auroral  field  line, 
the  dominant  collisional  process  in  the  ionosphere,  where  the  neutral- 
particle  density  greatly  exceeds  the  plasma  density,  is  ion-neutral  , 
charge  exchange.  The  collision  time  for  this  process  is  given  by 
t = l/ngov,  where  n^  is  the  neutral-particle  density,  a is  the  ion- 
neutral  charge- exchange  cross  section,  and  v is  the  magnitude  of  the 
relative  impact  velocity.  For  hot  magnetospheric  plasma  one  identifies 
V with  the  speed  of  the  charged  particle  and  requires  that  the  particle 
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be  relatively  free  of  collisions  over  its  full  bounce  period  2 tt  / Q-  ~ 5LR_/v 
Thus,  the  level  corresponding  to  a = X is  to  be  determined  from  the 
condition  that 


~ 2tt/ [ 5 nQ(f)  (rLRg  ] ~ 1 . 


(2) 


Using  an  average  CIRA  (1972)  model  of  the  neutral  atmosphere  and  a 

- 1 5 2 

representative  value  ( « 10  cm  ) for  the  proton-neutral  charge- 
exchange  cross  section  a,  we  determine  that  the  location  8=1  cor- 
responds to  an  altitude  ~2000  km.  The  corresponding  thermal  plasma 
density  at  s = / in  the  auroral  region  is  extremely  variable,  since 
sharp  ionospheric  troughs  are  not  limited  to  the  plasmapause  region 
(Hoffman  et  al.  , 1974).  At  high  invariant  latitudes  ( > 65°  ),  the 

trough  ionospheric  densities  may  be  as  low  as  10  -10  cm"  in  the 

4 - 3 

nightside  ionosphere  and  as  high  as  10  cm  in  the  dayside  ionosphere 
at  the  1400  km  altitude  (Hoffman  et  al.  , 1974).  While  the  criterion  (2) 
assures  us  that  magnetospheric  ions  will  be  relatively  free  of  charge- 
exchange  collisions  on  their  bounce  time  scale  in  the  interval  0 < a s I , 
charge-exchange  collisions  near  the  foot  of  the  field  line  may  exert  suffi- 
cient frictional  force  on  H^  ions  of  ionospheric  origin  to  modify  their 
density  distribution  along  the  field  line.  We  acknowledge  that  most 
energetic  eletrons,  as  well  as  ions  with  energies  i 50  keV,  remain 
collisionless  to  much  lower  altitudes  than  2000  km.  It  is  a shortcoming 
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of  our  model  that  all  particles  going  beyond  s = Z are  considered  lost. 
This  shortcoming  could  be  remedied  only  by  allowing  I to  depend  upon 
particle  species  and  energy,  but  such  a remedy  would  render  the  present 
model  entirely  intractable. 

Having  defined  the  interval  0 < s 5 for  which  our  model  applies, 
we  next  consider  the  consequences  of  Liouville's  theorem.  Since  the 
plasma  is  considered  to  be  collisionless  throughout  the  interval  OsssX.  , 
the  time- independent  distribution  function  for  either  species  must  depend 
only  on  the  constants  of  the  motion,  namely 

W = (n^q/2)(v„f  + vj)  + qlelV^  (3) 


and 


y- 


m V ^/2B 
q Is'  s 


(4) 


where  m^  is  the  particle  mass  and  q(=  ±1)  is  the  particle  charge  in  units 
of  the  electronic  charge  |ej  . Liouville's  theorem  asserts  that 


^±<"118' 


i 8 


f,  ( V , 

± IIS 


is' 


(5) 


if  the  points  (v  » v 5 s ) and  (v  ,,  v,  s'  ) are  connected  by  a 

Ro  XS  HS'  Xo 

dynamical  trajectory  in  phase  space.  The  conservation  laws  expressed 
by  (3)  and  (4)  imply  that 
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+ fl  - )]v  + (2qle!/m  )(V  - V ) 

B a i s ass 


for  any  such  pair  of  points  on  the  same  dynamical  trajectory  in  phase 
space.  Thus,  we  say  that  the  point  (v  , v ; s ) is  accessible  from  the 

MS  X 8 

source  point,  defined  as  (v  * , v * , s*  ),  if  v is  real  (i.  e.  , if 

•IS  X S il  S 

2 :k 

>0  ) for  all  s'  between  s and  s^  . An  obviously  necessary  con- 
dition for  particle  accessibility  to  s from  a source  point  s*  is  that 

V||^i  , as  given  by  (6),  be  positive  or  zero.  The  sufficiency  of  this 

2 

condition,  namely  v * 2 0,  depends  on  the  functional  form  of  V . 

II 8 3 

It  can  be  shown  (see  Appendix  A)  that  the  condition  v ^ > 0 is  sufficient 

It  s 

f ... 

for  accessibility  to  any  source  point  s for  both  species  of  particle 
( q = ±1 ) if  dV  , /dB  , >0  and 

S 3 


(d^V  ,/dB  ^,)  < 0 

s s 


for  all  points  s'  between  s and  s . 


In  our  model. 


we  use  the  criterion  v * 2 0 to  determine  particle  accessibility  to 

It  s 

their  respective  source  points  s ; thus,  we  require,  in  accordance  with 
(7),  that  V increase  monotonically  with  B and  appear  concave  down- 
ward  when  plotted  as  a function  of  B^  throughout  the  entire  interval 
Bq  S B^  S B^  . Since  the  functional  form  of  becomes  evident  only 
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upon  solution  of  (1),  we  have  to  assume  that  the  above  requirements  on 
2 2 

dV  /dB  and  d V /dB  are  satisfied  in  specifying  analytical  forms 

8 S S S 

for  f . (v  , V : V ) in  our  model.  Only  after  "evaluation"  of  the 
* II S 1 s s 

density  moment  and  "solution"  of  (1)  for  are  the  derivatives  dV^/dB^ 

2 2 

and  d V /dB  available  to  be  tested  for  sign.  These  conditions  thus  turn 
s s 

out  to  be  very  restrictive  constraints  (applied  a posteriori)  on  the  accep- 
tability of  "solutions"  to  (1).  In  theory,  any  "solution"  of  (1)  not  satisfy- 
ing these  constraints  has  to  be  discarded,  since  it  would  have  been  based 
on  a false  mapping  of  f±  in  (1).  In  practice,  however,  a slight  violation 
of  (7)  comparable  to  the  limit  of  numerical  resolution  in  computing  the 
density  moments  may  have  to  be  tolerated. 

A proper  consideration  of  constraints  upon  , based  on  accessi- 
bility of  a particle  to  its  source  point  in  phase  space,  should  include  not 
only  the  electrostatic  potential,  as  in  (7),  but  also  the  gravitational 
potential.  We  have  chosen  to  neglect  the  effects  of  gravity  in  (7)  because 
we  seek  solutions  for  which  le|V^  keV.  For  ions,  the  gravita- 

tional escape  energy  from  the  2000  km  altitude  is  ~8.2eV;  therefore,  we 
should  have  included  the  gravitational  potential  in  the  consideration  of 
phase-space  accessibility  of  cold  ions  near  the  foot  of  the  field  line  where 
|e|(  V,  - V ) 5 10  eV.  We  have  properly  included  the  gravitational 

potential  in  the  treatment  of  the  dynamics  of  cold  ions,  as  was  done  in 
the  works  of  Eviatar  et  al.  (1964)  and  Lemaire  and  Scherer  (1973).  How- 
ever, the  enforcement  of  (7)  as  a constraint  upon  near  the  foot 
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of  the  field  line  mey  not  be  aufficient,  i.  e.  , may  constitute  an  erroneous 
allocation  of  points  (v^,  v^;  s ) in  phase  space  with  respect  to  accessibility 
of  cold  ions  from  their  source  point  8*  = X.  We  have  b een  able  to  show 
that,  even  if  gravity  is  taken  into  consideration,  cold  ions  satisfying 
V i >0  in  (6)  are  accessible  to  their  source  at  the  foot  of  the  field  line 

H B 

if  V satisfies  (7)  and  if  dV  /dB  is  sufficiently  large  near 
the  foot  of  the  field  line  (see  Appendix  A). 

In  the  foregoing,  we  have  concentrated  on  discussing  the  properties 
of  imposed  by  the  requirement  that  the  particles  in  phase  space  be 

accessible  along  dynamical  trajectories  to  their  sources,  which  are 
assumed  to  be  either  at  the  equator  (s  = 0 ) or  at  the  foot  of  the  field 
line  (s  - X).  Among  the  possible  particle  populations  supporting  , 
there  may  be  populations  which  are  not  accessible  to  any  source  points 
in  the  interval  Os  s (Lemaire  and  Scherer,  1971  a;  b ).  Such  particles 
execute  dynamical  trajectories  with  two  turning  points  in  the  interval 
0 < 8 < X , and  are  therefore  trapped  by  the  electrostatic  potential  on  the 
one  hand  and  by  the  effective  magnetic  mirroring  potential  on  the  other. 

For  a potential  V ( 0 < V S V,  ) satisfying  the  constraint  imposed  by 
(7),  only  electrons  may  execute  such  trapped  trajectories.  We  assume 

that  the  phase- space  trajectories  available  for  such  trapped  electrons 

2 2 

are  defined  by  the  requirements  that  v^q  <0  and  v^^  < 0 . 

Details  on  the  construction  of  the  model  particle  distribution 
functions  for  the  various  populations  are  given  in  the  next  section. 
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PARTICLE  POPULATIONS 


In  the  present  model,  the  following  plasma- constituent  distributions 

are  assumed  to  be  present  in  the  collisionless  interval  ( 0 < s ^ -t  ) of  the 

auroral  flux  tube:  hot  anisotrop  c magnetospheric  plasma,  f.  . (v  v V ); 

Mq  lls  ’ is  ’ s 

cold  thermal  ionospheric  plasma,  f ( v.  ,v  ;V  );  two  populations  of 

X O S X s s 

backscattered  electrons,  fg.  ( ) with  ( i = 1,2);  and  trapped 

electrons,  f^  ( ^s  magnetospheric  distributions  fj^q  are 

considered  to  arise  from  an  equatorial  source  (s*  = 0 ) while  the  iono- 
spheric distribution  f,  and  the  backscattered  distributions  f„.  are  con- 

tq  Si 

sidered  to  be  injected  at  s = ^ . The  trapped  electrons  f,j,  are  considered 
to  have  been  scattered  into  their  trapped  trajectories  and  do  not  require  an 
accessible  source. 

2 

Particles  belonging  to  fj^q  are  classified  as  precipitating  ( V||^  >0  ) 

2 

or  mirroring  ( 0 ) according  to  their  respective  values  of 


II  8 


+ [1  - (B^/B^)]vj^^  + (2q|el/m  )(V  - V ), 


(8) 


which  determines  accessiblity  to  s = 1.  Since  V||  is  positive  for  down- 
going particles  in  our  convention,  and  since  all  magnetospheric  particles 
going  beyond  b = t are  considered  to  be  absorbed,  we  require  that  f.,  = 0 

for  trajectories  that  obey  both  the  condition  V||  <0  and  the  condition  v^^  > 0 . 

Further,  since  the  source  of  the  fj^^  is  located  at  s = 0,  we  require  that 


t. 
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all  trajectories  of  be  accessible  to  s = 0 by  the  criterion  discussed 

in  the  previous  section,  i.  e.  , that  ^0.  Noting  that  magnetospheric 

particles  may  be  anisotropic,  we  postulate  an  otherwise  bi- Maxwellian 

form  for  f.,  at  the  equator,  so  that 
M± 


f (v  V V ) 
Mq'  II s'  is'  s 


X exp  ) - (m^v,|^ /ZhT,,^)  - (q  I e I /k 
- - (Bq/B^)]  - (-qvJ/2xT,^)(Bo/B^)}  , (9) 


where  x is  the  Boltzmann  constant  and  is  a normalization  constant 

defined  in  such  a way  that  the  equatorial  number  density  is  given  by 


/ + C0  CO 

/ '^1  *^''^1  . V ; 0^ 

II  y i i Mq'  II  i 


The  unit  step  function  9(x)  is  defined  as  + 1 for  x 20  and  vanishes  for 


X <■  0 . 

For  ionospheric  electrons,  we  require  accessibility  from  the  source 

2 2 

at  8 = / , i.  e.  , we  require  that  V||^  >0  . The  condition  v^q  ^ 0 dis- 

2 

tingushes  between  electrons  that  mirror  before  reaching  the  equator  ( v^q  < 0 ) 

2 

and  those  that  cross  the  equator  { Vjjq  >0  ) to  be  lost  in  the  conjugate  iono- 
sphere. However,  the  "loss  cone"  is  assumed  to  be  completely  filled  in 
either  case,  and  the  distribution  function  is  not  complicated  by  this  dis- 
tinction. All  cold  electrons  that  enter  the  oaosphere  are  assumed  to  be 


k 


replenished  immediately  and  the  ionosphere  is  considered  completely 
symmetrical  between  north  and  south.  Thus,  we  postulate  an  iono- 
spheric electron  distribution  function,  which  is  taken  to  be  an  isotropic 


Maxwellian  distribution  of  temperature  > at  s = x 


fr  (v  , V : V ) 
1-  IIS  is  s 


X exp  I - (m_v^/2KTj_)  - (q  [ ej  ) | , (11) 

where  Nj . is  the  cold  ionospheric  electron  density  at  s = i . 

For  cold  ions,  the  effect  of  gravity  must  be  taken  into  account. 

Thus,  the  mapping  relationship  between  the  source  point  ( v , v I ) 

Ik  X XX 

and  the  point  ( 'll ^ s . s )<  analogous  to  (6),  is 


♦ [•  - t (2lel/m  l(V^  . V,) 


+ 2GMg[(l/r^)  - (1/rJ  ] , (12) 


where  G is  the  g ravitationcil  constant,  M^,  is  the  mass  of  the  earth, 

and  r is  the  radial  distance  between  the  center  of  the  earth  and  the 
s 

point  s , The  quantity  GM^m^^  ^ approximately  10.53  eV; 

therefore,  the  influence  of  gravity  upon  the  distributions  of  cold  ion- 
density  and  parallel  electrostatic  potential  is  important  near  the  foot 
of  the  field  line,  where  jej(  ) is  of  comparable  magnitude. 
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The  dominant  ionic  species  in  the  topside  ionosphere  are  O and  H 
The  mass  difference  between  the  two  species  influences  their  relative 
distributions  along  the  field  line  and  creates  an  additional  charge- separation 
electric  field,  as  in  the  polar  wind  (Banks  and  Holzer,  1968),  although  we 
are  dealing  with  closed  field  lines.  While  such  an  effect  may  be  important 
in  connection  with  the  distribution  of  and  ions  at  the  foot  of  the 

field  line,  we  note  again  that  such  a charge- separation  electric  potential 
IS  orders  of  magnitude  smaller  than  , which  is  of  the  order  of  kilo- 
volts. Thus,  for  the  sake  of  simplicity  in  exposing  what  are  possibly  the 
major  factors  in  the  maintainence  of  , we  assume  that  the  ionosphere 
consists  of  H ^ ions  only,  A treatment  of  the  mutual  consistency  between 
the  large  auroral  parallel  electrostatic  potential  and  the  distributions  of 
and  H ions  in  the  trough  region  of  the  ionosphere  will  be  given  in 
subsequent  work  on  the  formation  of  the  ionospheric  trough.  For  the 
present  paper,  we  use  the  cold  ion  distribution 


I + II  s is  3 


X ,xp  j-  v./ZxTitl  - . V,) 

- (GM^m^/KTj^)[(l/r^)  - (1/r^)  ] I , (13a) 


which  is  the  mapping  of  a simple  Maxwellian  distribution  in  accordance 
with  (12).  This  distribution  is  similar  to  that  postulated  by  Lemaire  and 


1 


► 


i 

i 

; 

I 

I 


Scherer  (1973):  we  assume  that  upward- moving  ions  are  lost  as  they 
cross  the  equator  so  that  there  are  no  downward-moving  ions  in  the 
interval  0 < s < i . 

For  the  sake  of  completeness,  we  have  also  experimented  with 
forms  for  which  the  ion  stream  has  speed  u ( > 0 ) at  s = / , as  in  the 
polar  wind: 


4+^"iis'"is=^s)  = v^i^)6(v  + u) 


X exp(- 


(13b) 


and 


(m^/4-n-K  ^1  + exp(-  m^  ® 

X [6(v^^^  + - u)e(v^^^)  j . 


(13c) 


The  form  (13b)  resembles  (13a)  and  that  used  by  Lemaire  and  Scherer 
(1973)  in  that  upward- moving  ions  are  considered  to  be  lost  as  they  cross 
the  equator.  The  form  given  by  (13c)  resembles  that  postulated  by  Schulz 
and  Koons  (1972)  in  that  it  consists  of  two  oppositely  directed  streams  at 
all  points  along  the  field  line.  The  mapping  of  (13b)  and  (13c)  is  in  accor- 
dance with  (12):  therefore,  we  have  to  make  sure  that  the  ion  stream  will 
be  accessible  to  s = t.  We  have  shown  that  ions  are  not  trapped  (so  as 
to  be  inaccessible  to  both  s = 0 and  s = t)  by  the  combined  electrostatic 
and  gravitational  potential  if  V varies  sufficiently  with  B near  the  foot 
of  the  field  line.  In  a proper  treatment,  utilizing  the  form  (13a)  for  the 
and  species  separately,  the  artificial  modeling  of  an  upward 

flow  at  the  foot  of  the  field  line  is  unnecessary. 
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In  the  construction  of  ion  distributions  above,  we  have  ignored 
the  important  effect  of  the  magnetosphe r ic  convection  electric  field  upon 
the  cold  ions  in  the  auroral  flux  tube.  Under  the  action  of  the  perpen- 
dicular convection  electric  field,  the  particles  begin  to  drift  across  B . 
Electron  drift  is  negligible  on  the  bounce  time  scale,  but  cold- ion  drift 
IS  not.  The  net  result  is  a loss  of  cold  ions  from  the  flux  tube.  Whereas 
Lemaire  and  Scherer  (1  973)  postulated  the  loss  process  to  have  occurred 
entirely  at  the  equator,  we  have  attempted  to  model  the  depletion  of  ions 
along  the  flux  tube  by  assuming  that  the  convection  electric  field  causes 
a constant  probability  of  cold- ion  depletion  per  unit  field-line  length. 
Thus,  the  distributions  (13)  are  to  be  multiplied  by  a depletion  factor 


where  t,  the  time  constant  for  the  process,  need  not  be  very  long  in 
comparison  with  the  cold- ion  bounce  period.  In  the  present  model,  we 
assume  t to  be  an  adjustable  parameter,  and  x is  generally  found  to  be 
of  the  order  of  20  minutes.  Since  the  inclusion  of  (14)  in  the  distribution 
function  (13)  of  cold  ions  renders  the  density  moment  integrations  entirely 
intractable,  we  assume  v„  , in  (14)  to  be  given  by  the  parallel  velocity 
of  a cold  ion  stream  free  of  the  influence  of  the  magnetic  field,  i.  e.  , by 
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-t-  (2  |e  I /m  )(V  - V ,) 

+ r s 


2GM  [d/r  ) - (1/r.)]  } 


1/2 


(15 


where  u,  assumed  to  be  approximately  one  half  of  the  cold- ion  thermal 
speed,  is  the  upward  flow  speed  at  the  foot  of  the  field  line.  With  the  use 
of  a trial  model  of  , , the  depletion  factor  D is  thus  made  indepen- 

dent of  the  phase- space  integrations  and  the  cold- ion  density  is  given  by 
the  product  of  the  density  moment  obtained  from  (13a)  and  the  depletion 
factor  D. 

Backscattered  electrons  are  assumed  to  be  produced  at  altitudes 

below  s = I in  a manner  that  yields  isotropic  distributions  at  s = I , 

but  3-1  is  regarded  as  the  source  of  backscattered  electrons  in  the 

=1=  2 

sense  that  accessibility  to  s = s requires  > 0 . As  in  the  model 

of  Evans  (1974),  which  is  based  on  the  abstracted  results  of  a detailed 
calculation  by  Banks  et  al,  (1974),  we  distinguish  two  populations  of 
backscattered  electrons:  the  primary  backscattered  electron  distribu- 
tion fc  , ( V ,y  ,V  ) and  the  secondary  distribution  fr, ( y,  , v,  ,V  ). 
Both  distributions  of  backscattered  electrons  are  related  to  the  precipi- 
tating magnetospheric  electron  flux  per  unit  kinetic  energy  s = I , 

namely 
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Zrr  f (2c, /m_^)  cos  o,  ^ ^ - V,  ) 6 (V||  ,^)  d(co9  o,  ) 

•'A  ’ ^ 


2TT(e,/m;)Cj^_  e(c,  - e V,  ) 


X exp[  - + ( |e  /kT||  _)  ] 


(16) 


We  estimate  the  relationship  between  the  upward  primary  back- 
scatter  flux  dJgj|/de^  and  the  total  precipitating  flux  dJ^/  de^  of 
energetic  electrons  by  means  of  a transfer  function  S (c,  , c,’  ) in 
the  integral  equation 


dj3„/d. 


S(e, , e,')  (d  /de,' ) 


(17) 


The  total  precipitating  flux  dJ.  /d  e/  in  (17)  is  given  by  the  sum  of 
the  magnetospheric  precipitating  flux,  (16),  and  the  precipi- 

tating flux  of  the  primary  backscattered  electrons,  dJgj^/de^’  . Kine- 
matical  considerations  require  that  the  discussion 

on  thermal  electrons  applies  equally  well  to  backscattered  electrons: 
whatever  goes  up  must  come  down.  For  the  primary-backscatter 
transfer  function  S (c,  , c^'),  we  choose  a form  given  by 

S(Cj,,5;)  = f €,/(e;)^]e(c;  - e^)  . (18) 
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This  represents  an  analytically  convenient  fit  to  the  transfer  function 
deduced  by  Evans  (1974)  from  a detailed  calculation  by  Banks  et  al. 
(1974).  The  normalization  of  (18)  yields  an  integrated  probability  of 
0.5  for  an  incident  electron  to  be  backscattered  with  energy  be- 
tween 0 and  . This  agrees  well  with  the  estimate  by  Evans 
(1974)  that  46%  of  the  incident  electron  flux  is  backscattered.  Thus  , 
we  obtain  the  relationship 


dJsp/dej^ 


= f fVK'‘ 


]f(dJ 


(dJsi,/dc 


dc. 


(19) 


where  we  have  expressed  J,  = J.,,  + Jc- 1 The  solution 

i M*  b t*  Mi  b IT 

of  (19)  for  in  terms  of  effected  by  multiplication 

of  both  sides  by  1 / and  differentiation  with  respect  to  . These 
steps  yield 


1 d 


d J 


ST 


d€. 


1 d J 


~T 


Mi 


de. 


The  solution  of  (20)  is  given  by 


00 

dJsi,/d€,  = J (dCj;/£i)(dJ„,/d€^) 


(20) 
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« (2TTKT|^_/m_^)Cj^_  |e(!e(  Vj^  - 

+ 0(Gj^  - |e|  Vjj^)  exp[(|e|  - ej^)/KT^^_]}  (21) 


upon  insertion  of  (16)  above.  The  integration  of  (16)  and  (21)  with  respect 
to  e,  reveals  a total  precipitation  flux 

= •^Mi  •^Sli  ” ^■^Mi  ^-^Sli 

« 4TT(HT|i_/m_)^Cj^_  f 1 + (|e|  Vj^/hT||_)]  . (22) 

I hus,  we  find  that  half  of  the  precipitating  energetic  electron  flux  is 
magnetosphe ric,  while  the  other  half  consists  of  electrons  previously 
backscattered. 

The  secondary  bac^scatter  is  characterized  (Banks  et  al.  , 1974; 
Evans,  1974)  by  a "universal"  spectrum  having  an  intensity  proportional 
to  the  total  flux  of  incident  energetic  electrons.  The  secondary  flux 
IS  negligible  (S  1%  of  maximum)  at  energies  e.  ~ 200  eV,  and  so  we 
obtain  a fit  in  simple  analytic  form  of  the  lower- energy  portion  of  the 
"universal"  secondary  backscatter  spectrum  given  by  Banks  et  al.  (1974): 

dJg2f/dej|^  **  ^ ^ 2^ 
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where  » 3.  0,  s:  15  eV,  and  Jj  is  given  by  (22).  In  the  application 

of  the  present  model  to  various  observations  to  be  considered  in  the  next 

section,  we  found  it  necessary  to  vary  somewhat  in  order  to  fit  the 

low-energy  portions  of  observed  electron  spectra;  therefore,  will 

be  considered  an  adjustable  parameter  of  the  model. 

It  remains  necessary  to  convert  the  upward  differential  fluxes 

given  by  (21)  and  (23)  into  distribution  functions  f„.  ( y,  ,v,  ;V  ), 

01  IIS  is  s 

where  i = 1,  2.  We  do  this  by  assuming  that  the  backscatter  distribu- 
tions are  isotropic  (cf.  Banks  et  al.  , 1974;  Evans,  1974)  at  s = l . One 
then  obtains  the  relationship 

dJsi,/dej^  ~ (2TT€jj^/m_^ ) fg^(vj|^  ;Vj^)  , (24) 


from  which  we  derive  the  following  backscattered  electron  distribution 
functions: 
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(2KT,,_/m_)Cj^_  e(v„^)|e(2|e(  - m_v^) 

+ e(m  v^  - 2 I e I V ) exp  [ ( 2 |e  I V - m v^)/2kT 

- a ' JL  a — a li 

T [ Vg^  + (2)e)  /m_)(V^  - V^)  ] (25) 


and 
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[1  t (ie|Vj/HT,  )]e(v,^) 

T [1  f (lel/e^HVj^  - V^)  + {m_v^/2r,2)f 

(26) 

For  electrons  oi  trapped  trajectories,  we  construct  their  distri- 
bution f_  ( V ,v.  ,V  ) by  as8umiii,i  it  to  be  a Maxwellian  distribution  in 
equiliorium  with  the  "source”,  from  which  the  trapped  electrons  had  been 
scattered  onto  their  present  trajectories  (Lemaire  and  Scherer,  1971b). 

In  our  model,  there  are  two  such  possible  "sources",  namely  those  at 
s = 0 and  at  s = £.  Thus,  for  each  "source"  j , we  construct  a trapped 
electron  distribution 

,27, 

X exp  j - (m  v^/ZkT)  4 (|e|/HT,)(V  - V *)  }, 

I - 3 J J 3 3j 

where  b.  is  the  location  of  the  "source"  and  T.  is  an  adjustable  tempera- 
ture  for  the  trapped  population.  For  j = 1 , we  identify  f,j,  ^ with  the 
distribution  scattered  into  trapped  trajectories  from  fj^  , and  so  we  take 
Nj  ( s J ) = Nj^  . For  j = 2 we  identify  f,p  ^ with  the  distribution  of 
electrons  that  have  been  scattered  into  trapped  trajectories  from  distri- 
butions whose  sources  are  at  a = I . There  are  several  such  "source" 
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distributions,  but  calculations  reveal  that  the  density  moments  of  low- 
temperature  trapped  particles  originally  injected  at  a - L are  negligible 
due  to  phase-space  limitations  if  V ~1  kilovolt.  Therefore,  we  identify 
with  electrons  scattered  into  trapped  trajectories  from  fg  ^ . Since 
the  primary  backscattered  electron  distribution  has  temperature  , 

we  assume  both  trapped  distributions  to  have  the  same  temperature 
T = T I , forming  a single  distribution 


(m_/2TiKT)^^^N.j,  e(-V|,^)  e(-vj) 

X exp  [ - (m  v^^/2hT)  + ( e V /^T)]  . (28) 


where 

" ^M-  ^ Ngjexp(  - |e|  V^/hT)  , (29) 

Ng  j being  the  density  derived  from  fg  ^ at  a - I . The  existence  of 
trapped  electrons  had  been  assumed  in  previous  exospheric  models  of  the 
solar  wind  (Lemaire  and  Scherer,  1971b),  but  not  in  the  magnetosphere 
(Eviatar  et  al,  , 1964;  Lemaire  and  Scherer,  1973).  Observational  evi- 
dence for  or  against  their  existence  would  be  welcome.  Perhaps  one 
may  argue  that  certain  regions  of  phase  space  would  be  empty  if  trapped 
electrons  do  not  exist  (see  Figure  lb).  Evidence  suggesting  the  filling 
of  all  regions  of  phase  space  in  auroral  electron-precipitation  events  is 
given  by  Kaufmann  et  al.  (1976). 


Fif?ure  1 summarizes  the  regions  of  phase  space  occupied  by  the 
various  particle  populations  discussed  above.  The  left  panel  indicates 
veloc ity- space  demarcations  that  characterize  the  various  portions  of 
the  magnetospheric  and  ionospheric  populations  of  ions.  The  right 
panel  illustrates  the  corresponding  velocity-space  demarcations  for 
electrons. 
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Figure  1.  The  regions  of  phase  space  occupied  by  the  various 
particle  populations  are  shown  schematically  for  (a)  ions  and  for  (b) 
electrons.  The  solid  curves  are  demarcations  in  phase  space  for  the 
various  populations  of  particles  identified  by  labels  in  circles  or  in 
cartouches:  M for  particles  of  magnetosphe ric  origin,  1 for  particles 

of  ionospheric  origin,  S for  backscattered  electrons,  and  T for  electrons 
trapped  between  mirror  points  that  both  lie  on  the  same  half  of  the  field 
line.  The  ionic  label  I in  parentheses  refers  to  the  population  given  by 
(13c),  as  distinguished  from  those  given  in  (13a)  and  (13b).  The  elec- 
tronic labels  1 in  parentheses  indicate  phase  space  regions  occupied  by 
electrons  in  the  extreme  tail  of  the  cold  ionospheric  energy  distribution. 
The  dashed  diagonal  lines  are  the  asymptotes  v^  = ± [ (B^  /B)  - 1 ] ^ v^^ 

of  the  hyperbolic  phase  space  demarcations  between  v^^^  > 0 and  v^^^  < 0. 


THE  MODEL 


In  the  previous  two  sections  we  have  summarized  the  manner  in 
which  charge  neutrality  is  to  be  enforced  on  the  constructed  particle 
distributions  in  an  auroral  flux  tube.  The  general  features  of  the  model 
dealt  with  in  the  previous  sections  are  quite  akin  to  those  of  the  exo- 
sphere and  plasma- sheet  models  of  Lemaire  and  Scherer  (1971b,  1973, 
1977),  although  considerations  such  as  accessibility  of  sources,  anisotropy 
of  magnetospher ic  particles,  and  the  inclusion  of  backscattered  electrons 
characterize  the  greater  generality  of  the  present  formulation.  However, 
it  is  mainly  the  area  of  model  usage  which  distinguishes  the  present  model 
from  previous  exospheric  models,  since  the  major  purpose  of  the  present 
model  is  to  explore  the  factors  that  may  contribute  to  the  maintainence 
of  a large  auroral  potential  difference  between  the  equator  and  the  iono- 
sphere. Therefore,  we  apply  our  model  to  specific  electron-precipitation 
events  for  which  distinctly  beam-like  energetic  electrons  are  observed. 

A second  feature  involving  the  usage  of  our  model  has  to  do  with 
the  geometry  of  the  assumed  magnetic  field.  In  the  construction  of 
particle  distributions  in  the  previous  section,  we  have  tacitly  assumed 
that  the  model  is  to  be  applied  to  particle  and  electric-field  distributions 
in  a closed  dipolar  flux  tube.  While  we  shall  continue  to  assume  so  in 
the  rest  of  this  work,  we  wish  to  note  that  the  physical  processes  dis- 
cussed in  this  work  are  also  applicable  to  regions  of  open  field  lines. 

For  regions  of  open  field  lines,  it  would  be  necessary  to  reconstruct 
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the  particle  distributions  discussed  in  the  previous  section.  Although 
magnetospheric  particles  in  such  open-field-line  regions  are  likely  to 
be  isotropic,  a large  electrostatic  potential  difference  may  still  be 
maintained  by  virtue  of  electron  backscatter.  Applications  of  our 
model  of  particle  and  electric  field  distributions  for  open  field-line 
regions  will  be  considered  in  a subsequent  work. 

In  typical  observations,  either  the  electron  flux  distribution  (at  a 
given  pitch  angle)  dJ_/dedil  or  the  electron  energy  flux(at  a given 
pitch  angle)  «:dJ_/dedi2  is  observed  at  an  altitude  below  2000  km.  As 
a first  step  in  the  usage  of  our  model,  we  require  that  the  theoretical 
flux  distribution  at  the  foot  of  the  field  line,  i ‘f., 

(dJ_/dc^dU)  = + ^S1  ® 

where  a is  the  pitch  angle  and  dS2  the  corresponding  element  of  solid 
angle,  be  fitted  to  the  observed  flux  distribution.  Except  for  events  in 
which  diffusive  energy  gain  is  evident  (Evans,  1975),  good  fits  are  usually 
obtained.  The  parameters  , ( Nj^  /T^  ) , and  are  strin- 

gently constrained  by  the  fitting.  Specific  examples  of  fits  to  observed 
particle-flux  and  observed  energy-flux  distributions  are  shown  in  Figures 
2 and  3,  respectively.  Detailed  discussions  of  the  cases  will  be  given  in 
the  next  section.  In  the  usual  situation,  only  the  electron-flux  distribution 


-40- 


ENERGY,  keV 


Figure  3.  Electron  energy  fluxes  at  0°  and  at  180°  pitch  angles 
observed  in  an  "inverted- V structure  (Mizera  et  al.  , 1976)  are  shown 
as  functions  of  electron  energy.  The  difference  between  the  downgoing 
and  upgoing  fluxes  at  1.76  keV  indicates  the  presence  of. a "monoenergetic 
beam  of  precipitating  energetic  electrons. 


for  a specific  event  is  available  for  comparison  with  the  theoretical  flux; 
therefore,  we  assume  the  ion  parameters  to  be  variable  in  the  search  for 
a self-consistent  parallel  electrostatic  potential  . Should  energetic 
ion  distributions  be  available,  a self-consistency  check  of  our  model  can 
be  made  immediately  since  is  a common  parameter  for  both  electron 

and  ion  distributions.  In  the  usage  of  our  model,  the  field- aligned  current 

CD 

V dv  / [ f , (v  , V : V ) 

is  is  J ^ f'  II s'  is  3 

- CD 

IS  neither  made  to  vanish  nor  assigned  any  particular  numerical  value, 
as  is  done  in  Lemaire  and  Scherer  (1973).  If  is  measured  for  a 
specific  event  to  which  we  apply  our  model,  then  an  additional  consistency 
ckeck  can  be  made  on  the  model  since  is  again  a common  parameter 

in  (30)  and  (31).  In  general,  because  of  the  assumed  magnetospher ic  particle 
anisotropy,  and  because  of  the  loss  of  cold  ions  to  the  magnetospheric 
convection  pattern,  is  nonvanishing  in  our  model.  Further,  the  loss 
of  cold  ions  from  the  flux  tube  implies  that  a perpendicular  current 
also  exists.  Such  considerations  raise  the  question  of  how  to  complete 
the  electrical  circuit  in  the  auroral  region.  This  presumably  occurs  via 
other  field  lines  than  the  one  under  consideration  cind  presumably  involves 
perpendicular  electric  fields  which  co-exist  with  the  parallel  electric 
field.  It  is  not  our  intention  to  explore  here  the  mechanisms  by  which 


- f (v, 


V ;V  ) ] 
is  s J 


dv 


II  s 
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any  such  magnetospheric  circuits  are  closed.  In  subsequent  work,  such 
questions  will  be  addressed  in  connection  with  the  so-called  "inverted-V" 
structures. 

Having  obtained  the  parjuneters  ^ ^ "^2 

from  a comparison  of  (28)  with  data  for  a specific  event,  we  next  proceed 
to  implement  the  quas  i -neut  rality  condition  (1).  Since  we  have  assumed 
specific  forms  for  the  distribution  functions,  we  calculate  the  density 
moments  of  ^I±'  ^Si  ^ closed  form  ) as  functions  of 

’^s’  ^s’  the  particle  parameters  ( , Nj^,  T||^  , Tj^,  T,  A^, 

and  ' ).  Examples  of  the  required  velocity- space  integrations,  some  of  which 
are  quite  tedious,  are  given  in  Appendix  B.  Not  all  of  the  above  particle 
parameters  are  entirely  independent,  since  the  limits  of  (1)  as  s -•  0 and 
as  B ->  L constitute  two  constraints  among  the  above  particle  parameters 
and  the  potential  difference  . Thus,  only  the  parameters  » 

> Tj^  , T,  and  t may  be  considered  free  parameters  of  the  model, 
although  the  known  conditions  of  the  ionosphere  in  the  auroral  region  impose 
additional  bounds  on  Nj  ^and  . A summary  of  the  parameters  of  the 

model  IS  given  in  Table  1. 

Once  the  particle  parameters  are  specified,  subject  to  the  constraints 
discussed  above,  we  consider  (1),  the  quasi- neutrality  equation  in  integrated 

form,  as  a transcendental  equation  to  be  solved  for  V as  a function  of  B . 

8 s 
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Table  1.  Parameters  of  the  Model 


■A 


This  is  done  on  the  CDC  7600  computer,  since  a large  number  of  tran- 
scendental functions  are  involved.  The  "solution”  so  obtained  is, 
however,  not  necessarily  the  required  solution  of  the  model,  since  an 
acceptable  must  be  a monotonic  function  of  and  it  must  satisfy 

(7).  If  does  not  satisfy  these  tests,  the  particle  parameters  are 

adjusted  and  the  solution  process  is  repeated  until  a satisfactory  solution 
is  obtained. 
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RESULTS 


Since  the  primary  purpose  of  the  present  model  is  to  determine 
circumstances  under  which  a large  electrostatic  potential  difference 
may  he  maintained  along  an  auroral  field  line,  we  have  applied  the  model 
to  a number  of  specific  cases,  two  of  which  will  be  discussed  in  this 
section. 

Figure  2 shtiws  the  unidirectional  electron  flux  at  4S*’  pilch  ancle 
for  an  event  observed  in  the  dayside  cusp  region  on  11  January  l'*7‘S 
(Winningham  et  al.  , 1977).  The  electron-flux  data  have  been  compared 
with  the  model  (solid  curve)  of  Evans  (1974), and  the  ficure  is  taken  from 
Evans  (1975).  The  data  suggest  that  a monoenergetic  beam  of  particle^' 
at  an  energy  ~600  eV  may  be  present  in  the  electron  spectrum.  Figure 
3 shows  electron  energy- flux  data  at  two  pitch  angles  obtained  in  an 
"inverted-V"  structure  (Mjy.era  et  al.  , 1976).  This  example  is  partic- 
ularly interesting  not  only  because  of  its  occurrence  in  the  auroral 
region  (invariant  latitude  69.  7*^’  ) but  also  because  it  clearly  shows  the 
beam-like  character  of  precipitating  high-energy  electrons,  as  is 
evidenced  by  the  great  difference  between  the  unidirectional  energy 
flux  at  0*^  and  that  at  180*^  pitch  angles.  The  parameters  which  deter- 
mine the  model  fits  shown  in  Figure  2 ( Case  W ) and  in  Figure  3 
( Case  M ) are  summarized  in  Table  2. 
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Table  2:  Parameters  of  the  Fit 


Parameter 


Case  W ( Figure  2 ) 


Case  M ( Figure  5 ) 


V 

i 


0.  59  kV 
5 cm 

0.  189  keV 
0.  775  keV 
4.  8 


1 . 76  kV 

0.  6 cm 

1 . 2^  keV 
4.  67  keV 
3.  2 


Altliough  the  parameters  shown  in  Table  2 are  not  uniquely  determined 
from  the  data,  the  comparison  between  the  present  model  parameters  and  those 
of  P'vans'  model  for  Case  W illustrates  some  differences  between  the  two 
models.  In  Evans'  model,  the  number  density  of  the  electron  beam  was  taken 

. 3 

to  be  0.  85  cm  and  the  temperature  was  taken  to  be  0,  143  keV.  While  his 
electron  temperature,  which  determines  the  slope  of  the  flux  distribution  at 
energies  above  '»‘iv  , , is  approximately  the  same  as  our  parallel  tempera- 
ture hT  , the  electron  beam  density  required  in  Evans'  model  is  consider- 
ably smaller  than  the  Nj^  of  our  model,  primarily  because  of  the  fact  that 
in  our  model  only  those  m agnetosphe ric  electrons  that  would  mirror  at  s >£ 
contribute  to  the  energetic  electron  flux.  Further,  because  of  differences 
in  the  pitch-angle  distributions  of  the  energetic  electrons  for  the  two  models, 
the  total  flux  (integrated  over  pitch  angle)  is  larger  in  Evans  model  than 
in  ours;  hence,  the  low-energy  electron  fluxes  for  the  two  models  are  also 
different.  For  Case  M (Figure  3),  the  required  electrostatic  potential 
difference  and  the  required  magnetospher ic  electron  temperatures  are  con- 
siderably higher  than  in  Case  W (Figure  2 ),  In  many  respects.  Case  M 
is  a more  stringent  test  of  the  model  since  the  electron  energy  fluxes  at 
two  pitch  angles  are  available.  In  particular,  the  low-energy  portion  of 
the  electron  spectrum  shown  in  Figure  3 requires  considerably  lower  secon- 
dary backscatter  than  in  Figure  2. 

With  the  use  of  the  parameters  determined  in  Table  2,  searches  for 
self-consistent  electrostatic  potentials  were  made  in  accordance  with 

the  formulation  given  in  the  previous  two  sections.  In  the  solution  process. 
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the  model  parameters  not  yet  determined  in  Table  2 are  varied  in  the  manner 
described  in  Table  1.  For  Case  M,  the  L value  of  the  field  line  is  8.35  . For 
Case  W,  the  observations  were  made  at  very  high  invariant  latitudes  (77.  5°  N 
where  the  rocket  was  possibly  on  open  field  lines  in  the  cusp  region.  In  this 
respect,  our  model  is  probably  not  fully  applicable,  since  we  assume  the 
magnetosphe r ic  particles  to  be  partially  trapped  by  a dipolar  field.  However, 
if  the  precipitating  magnetospheric  particles  were  to  have  the  anisotropies 
determined  in  the  fitting  of  Figure  2,  the  application  of  the  model  is  not 
substantially  changed  even  though  the  solution  V^,  as  a function  of  B/Bq  , 
may  not  be  an  accurate  description  of  the  "true"  potential  distribution. 

Since  the  requirement  of  the  total  potential  drop  remains  the  same 

whether  one  considers  open  or  closed  field  lines,  we  attempted  an  appli- 
cation of  our  model  to  Case  W by  assuming  an  L value  of  8.  54  for  this 
event.  For  this  reason,  our  results  for  Case  W are  rather  to  be  viewed 
as  an  illustration  of  our  model  them  as  reconstructions  of  actual  particle 
and  potential  distributions.  The  electrostatic  potential  solutions  of  (1) 
satisfying  the  accessibility  criterion  (7)  for  Cases  W and  M are  shown  in 
Figure  4 (Case  W)  and  Figure  5 (Case  M)  respectively.  In  Figure  4,  it 
IS  seen  that  the  potential  distribution  is  concave  downward  when 

plotted  against  the  magnetic- field  ratio  B/Bq.  In  Figure  5,  the  potential 
distribution  is  also  concave  downward  although  it  may  not  appear  so 

tiecause  the  magnetic-field  ratio  is  plotted  on  a logarithmic  scale.  For 
Case  W (Figure  4),  the  major  part  of  the  potential  drop  is  confined  to  the 
region  above  ~ 10000  km  altitude,  with  parallel  electric  field  intensities 


MAGNETIC  FIELD  RATIO,  B'Bq 


Figure  4.  The  self-consistent  electrostatic  potential  computed 
for  the  case  corresponding  to  the  observed  electron  flux  of  Figure  Z 
(Case  W)  is  shown  as  function  of  magnetic -field  ratio  and  altitude.  T1 
magnitude  of  the  parallel  electrostatic  field  is  also  shown. 
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Figure  5.  The  self-consistent  electrostatic  potential  computed 
for  the  case  corresponding  to  the  observed  electron  flux  of  Figure  3 
(Case  M)  is  shown  as  function  of  magnetic-field  ratio  and  altitude.  The 
curve  for  the  potential  V has  the  same  property  of  downward  concavity 
as  that  shown  in  Figure  4.  The  difference  in  appearance  between  this 
figure  and  Figure  4 is  due  to  the  logarithmic  scale  used  in  plotting 
the  magnetic- field  ratio  here.  The  magnitude  of  the  parallel  electro- 
static field  computed  for  this  case  is  considerably  larger  than  that  of 
Figure  4. 
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of  up  to  40  uV/m,  The  potential  difference  for  Case  M is  considerably 

larger  than  that  for  Case  W,  hence,  the  parallel  electric  field  intensity  in 
Figure  5 (up  to  140  uV/m  ) is  considerably  larger  than  that  in  Figure  4. 

Quite  interestingly,  the  region  of  major  potential  drop  (above  ~ 12  000  km) 
in  Case  M is  not  drastically  different  from  that  in  Case  W,  even  though 
the  magnitudes  of  the  potential  drops  are  quite  different.  The  self-con- 
sistent electrostatic  potential  distributions  for  both  cases  are  consistent 
with  Evans  model,  where  it  is  assumed  that  the  region  of  electron 
acceleration  by  parallel  electric  fields  is  well  above  the  region  of  electron 
flux  observation.  In  both  cases,  the  field  lines  below  10  000  km  altitude 
are  essentially  equipotentials.  This  common  feature  may  be  a weak  point 
in  the  model,  since  there  is  some  observational  evidence  that  particle 
acceleration  occurs  mainly  below  the  7000  km  altitude  (J.  F.  Fennell  and 
P,  F.  Mizera,  personal  communications,  1977).  To  counterbalance  the  point, 
however,  there  also  is  some  evidence  of  parallel  electric  fields  above 
7000  km  from  bar ium  - release  experiments  (Haerendel  et  al,  , 1976). 
Evidently,  the  observational  situation  is  probably  quite  variable.  On 
the  other  hand,  the  location  of  intense  parallel  electric  field  in  our  model 
may  be  somewhat  dependent  on  the  assumed  location  of  the  ionospheric 
boundary.  Presumably,  if  we  were  able  to  handle  the  proper  ionospheric 
boundary  condition,  which  depends  on  particle  species  and  particle  energy, 
the  distribution  of  might  be  somewhat  different.  Moreover,  the 

inclusion  of  an  intense  perpendicular  electrostatic  field  in  the  considera- 
tion of  charge  neutrality  would  also  affect  the  location  of  the  major  parallel 
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electrostatic  potential  drop.  Evidence  for  the  important  influence  of 
very  intense  perpendicular  electric  fields  upon  energetic  particle  accelera- 
tions in  the  auroral  region  has  recently  been  reported  (Mozer  et  al,  , 1977). 
Perpendicular  electrostatic  fields,  moreover,  may  play  a major  role  in 
defining  inverted- V structures  (e.  g.  , Swift,  1975). 

The  critical  model  parameters  required  for  the  solution  of  (1), 
satisfying  the  accessibility  criterion  (7),  are  shown  in  Table  3.  All 
other  model  parameters  are  either  fixed  or  determined  by  charge-neu- 
trality constraints  at  s - 0 and  s = x’  , in  accordance  with  the  summary 
in  Table  1.  Comparison  of  magnetospher ic  ion  temperatures  in  Table  3 
witi'  m agnetospher ic  electron  temperatures  in  Table  2,  for  both  cases, 
indicates  that  the  solutions  require  highly  anisotropic  electrons  ( T^  /T^ 

= 4.1  for  Case  W and  T /T  =3,8  for  Case  M ) and  isotropic  ions 
( T^^/T||^  = 1 ),  even  though  the  potential  differences  and  particle 

densities  in  the  two  cases  are  quite  different.  Although  this  peculiar 
feature,  common  to  both  cases,  may  be  a coincidence,  we  are  inclined 
toward  the  view  that  the  difference  between  the  pitch- angle  anisotropies  of 
electrons  and  ions  required  by  the  model  is  not  accidental.  For  field  lines  just 
poleward  of  the  plasmapause,  where  the  magnetospher  ic  particle  populations 
are  probably  of  ring-current  origin,  the  ion  population  is  susceptible  to 
cyclotron  instability  at  lower  pitch- angle  anisotropies  than  the  corre- 
sponding process  for  the  electrons  (Cornwall  et  al.  , 1971).  Thus,  one 
would  expect  the  equilibrium  population  to  be  much  more  nearly  isotropic 
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Talilu  3:  Parameters  Yielding  Solutions 


Parameter 

Case  W ( Figure  4 ) 

Case  M ( Figure  ‘ 

1 ZOO  cm  ^ 

00  cm"  ^ 

2.  0 keV 

«.  0 keV 

2.  0 keV 

K.  0 keV 

t 

1 7 min 

19  min 

for  protons  than  for  electrons.  For  field  lines  nearer  to  the  cuap,  where 
the  magnetic  field  is  very  weak,  the  pitch-angle  anisotropy  of  an  ion 
population  IS  more  likely  to  be  destroyed  by  scattering  off  magnetic- 
field  inhomogeneities  of  a given  scale  size  than  is  the  anisotropy  of  an 
electron  population.  Even  without  consideration  of  scattering  by  random 
magnetic-field  inhomogeneities,  the  motion  of  ions  at  a given  energy  is 
less  likely  to  be  adiabatic  than  that  of  electrons  at  the  same  energy. 

Thus,  cons  iderations  of  possible  wave-particle  interactions  and  of  non- 
adiabaticity  of  particle  motion  in  the  earth's  magnetic  field  seem  to 
favor  a tendency  for  auroral-zone  ions  to  be  more  nearly  isotropic  in 
pitch-angle  distribution  than  electrons.  Such  a tendency  is  consistent 
with  the  temperature  parcuneters  required  for  self-consistent  solutions 
in  our  model. 

An  important  feature  of  our  model  concerns  the  relationship  be- 
tween the  magnitude  of  the  parallel  potential  drop  V,  and  the  ionospheric 
dnesity  ^ required  to  obtain  a self-consistent  solution.  Comparison 
of  Case  W with  Case  M in  Tables  2 and  3 shows  that  V,  and  N,  , are 
inversely  related,  i.  e.  , to  maintain  V,  = 0,  59  kV  in  Case  W requires 

JL 

Nj  ^ - 1 200  cm"  ^ , whereas  to  maintain  1.76kV  inCaseM  requires 

. 3 

cm  . The  sensitivity  of  our  model  to  the  ionospheric  density 
Nj  ^ and  (in  particular)  the  above  inverse  relationship  are  expected, 
since  is  essentially  maintained  by  charge  separation  caused  by  the 
difference  between  magnetospheric  particle  pitch-angle  anisotropies  and 
Viy  the  inevitability  of  electron  backscatter.  The  existence  of  a parallel 
potential  distribution  implies  that  ions  from  the  ionosphere  will  be 
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accelerated  up  tlie  field  line  to  neutralize  the  electron  excess.  Thus, 
the  higher  the  ion  density  Nj  ^ , the  more  difficult  it  is  for  the  model  to 
maintain  a large  equilibrium  potential  difference  . Although  local- 
time dependence  of  the  ionospheric  density  is  not  introduced  into  the 
model,  the  values  of  Nj  ^ required  for  both  cases  in  Table  3 are  con- 
sistent with  the  local  time  of  observation  for  the  two  cases.  Since 
Case  \V  was  observed  in  the  dayside  cusp  region  and  Case  M was  observed 
at  a local  time  of  2024  hours,  the  anticipated  day-night  difference  in  iono- 
spheric densities  at  2000  km  is  consistent  with  the  order-of- magnitude 
difference  between  Nj  ^ for  the  two  cases.  It  should  be  noted  that  Nj  ^ 
was  solely  determined  by  the  requirement  of  obtaining  a self-consistent 
solution.  The  density  Nj.  ^ determined  for  Case  M is  consistent  with  the 
winter  nighttime  trough-region  ionosphere  (Hoffman  et  al.  , 1974).  The 
density  ^ determined  for  Case  W is  also  consistent  with  the  winter 
dayside  high- 1 atitude  ionosphere,  smce  the  event  was  observed  in  January 
1976  (Wmningham  et  al.  , 1977).  As  is  discussed  above,  the  maintenance  of 
a kilovolt-magnitude  para  lei  electrostatic  potential  requires  rather  low 
ionospheric  densities,  wh  ch  correspond  to  winter  conditions  of  the  polar 
F-regi'in.  Conversely,  only  very  small  parallel  potential  differences  may 

be  maintained  if  the  F-region  ionospheric  densities  ^ at  2000  km  were 
4-3 

~ 10  cm  , corresponding  to  summer  daytime  conditions.  These  con- 
siderations are  applicable  only  to  the  general  dependence  on  the  ionospheric 
boundary  conditions  of  our  model  since  the  dayside  north- south  asymmetry 
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of  the  ionosphere  is  neglected  in  our  consideration.  Quite  interestingly, 
statistical  analysis  of  OGO-4  data  indicated  that  field- al igned  electron 
precipitation  events  were  more  likely  to  occur  during  winter  (Berko  and 
Hoffman,  1974).  Since  our  model  requires  ion  densities  at  the  lower 
boundary  which  are  consistent  with  F- region  trough  densities  only,  and 
further,  since  energetic  O and  H ions  in  the  keV  range  (moving 
upward  along  auroral  field  lines)  have  been  reported  (Shelley  ct  al.  , 1976), 
we  may  perhaps  invoke  the  self-consistent  parallel  electric  field  as  a 
possible  mechanism  for  the  formation  of  the  ionospheric  trough.  A detailed 
consideration  of  such  a mechanism  will  be  made  in  a subsequent  work. 

The  density  distributions  m of  particles  of  various  species 
( j = M-  , M4  , SI,  S2,  T,  I-  , and  1+  ) for  Case  W are  shown  m 
Figures  6 and  7.  From  Figure  6,  it  is  seen  that  roughly  1/3  of  the 
magnetospheric  electrons  are  precipitating,  so  as  to  yield  a precipitating 
electron  beam  of  density  0.  9 cm"  , in  agreement  with  the  density  required 
in  the  model  of  Evans  (1975).  The  density  distribution  ng^  of  primary 
backscattered  electrons  shows  a peculiar  kink  near  the  magnetic-field 
ratio  B/Bq~  300.  The  formation  of  the  kink  is  quite  easily  understood 
if  one  recalls  that,  according  to  (21),  the  primary  backscattered  electron 
distribution  consists  of  a low-energy  component  ( < l‘'|V^  ) and  a high- 

energy  component  ( > K'lVj  ) . The  high-energy  component  of  primary 

backscattered  electrons  has  a density  distribution  similar  to  nj^_  near 
the  foot  of  the  field  line,  whereas  the  low-energy  component  has  a dis- 
tribution similar  to  the  n^_  of  secondary  backscattered  electrons. 
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Figure  6.  Distributions  of  the  various  energetic  particles  along 
the  field  line,  computed  for  Case  W,  are  shown  as  functions  of  magnetic- 
field  ratio.  The  subscripts  identify  the  various  particle  populations:  M± 

for  magnetosphe r ic  protons  and  electrons,  SI  for  primary  backscatte red 
electrons,  and  T for  electrons  trapped  between  mirror  points  that  both  lie 


on  the  same  half  of  the  field  line. 


shown  in  Figure  7.  The  sum  of  two  such  curves  produces  the  shape  of  ngj 
shown  in  Figure  6.  Comparison  of  magnetospheric  ion  density  in 

Figure  6 with  ionospheric  ion  density  n^  ^ in  Figure  7 in  the  region  just 
above  the  intense  parallel  electric  field  indicates  that  the  densities  of  the 
two  species  are  comparable  there.  Since  the  ionospheric  ions  in  this 
region  have  been  accelerated  by  the  electric  field,  we  would  expect  sub- 
stantial numbers  of  energetic  ions  in  the  region,  in  agreement  with 

Shelley  et  al.  (1976). 


J 


I 

» 
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CONCLUSION 


We  have  demonstrated  in  this  model  that  parallel  electrostatic  potential 
differences  of  the  order  of  one  kilovolt  may  be  self- consistently  maintained 
alonjj  an  auroral  field  line  by  the  magnetospheric  and  ionospheric  particle 
populations  in  the  flux  tube.  Our  model  attempts  to  abstract  the  characteristics 
of  magnetospheric  and  ionospheric  particle  distributions  and  to  incorporate 
the  very  important  process  of  electron  backscatter  in  a magnetic  field.  By 
application  of  the  principle  of  quasi-neutrality,  a self-consistent  parallel 
electrostatic  potential  distribution,  in  equilibrium  with  the  particle  distribu- 
tions along  the  auroral  field  line,  is  obtained.  Application  of  the  model  to 
two  typical  events  of  field- aligned  electron  precipitation  have  yielded  particle 
parameters  which  are  consistent  with  the  conditions  of  the  ionosphere  and 
the  magnetosphere  at  the  time  and  place  of  observation  of  the  events.  The 
difference  in  the  pitch-angle  anisotropies  of  magnetospheric  ions  and  electrons 
and  the  inclusion  of  backscattered  electrons  are  the  key  factors  which  self- 
consistently  support  a parallel  potential  difference  of  kilovolt  range.  Our 
model  suggests  that  such  large  parallel  potential  differences  between  the 
ionosphere  and  the  magnetospheric  equator  can  be  maintained  only  if  the  ion 
density  at  the  lower  boundary  (2000  km  altitude)  is  comparable  to  that  in  the 
F- region  trough.  We  are  inclined  to  view  the  generation  of  a large  parallel 
electrostatic  potential  drop  as  a self-consistent  mechanism  by  which  the 
F-region  ionospheric  trough  is  formed  by  extraction  of  cold  ions  from  the 
ionosphere.  These  cold  ions  are  accelerated  upward  into  the  magnetosphere 
where  their  energies  will  be  ~|p|  , i.  e.  , in  the  keV  range. 
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As  we  have  stated  in  the  introduction,  our  model  is  not  intended  to 

be  applicable  to  all  types  of  field- aligned  electron  precipitation  events  in 

which  beam-like  characteristics  are  evident.  In  particular,  our  model  is 

incapable  of  explaining  events  in  which  diffusive  energy  gain  is  evident. 

Further,  our  model  is  not  intended  to  exclude  the  possible  existence  of 

double  layers  or  electrostatic  shocks.  Indeed,  it  is  quite  easy  to  obtain 

in  our  model  a ' solution"  which  has  the  characteristics  of  a double  layer 

similar  to  that  obtained  by  Lemaire  and  Scherer  (1977).  However,  such  a 

potential  distribution  does  not  satisfy  our  accessibility  criterion  imposed 

by  (7).  We  do  not  know  how  to  treat  the  veloc ity- space  integration  for 

particles  that  are  accessible  to  the  source  point  s’‘'  only  under  conditions 

2 

more  stringent  than  the  simple  criterion  on  Vjig*  • Therefore,  we  are 
unable  to  generate  double  layer"  solutions  that  are  acceptable  within  the 
franiework  of  our  model. 

In  a simple  model  such  as  ours,  there  are  obviously  a number  of 
features  which  require  improvement.  Clearly,  our  assumption  of  an 
ionospheric  boundary  independent  of  particle  species  and  energy  is  not 
justified.  The  exclusion  of  perpendicular  electric  fields  in  our  considera- 
tion of  quas  1- neutrality  is  clearly  not  realistic,  since  auroral  events  are 
invariably  accompanied  by  intense  perpendicular  electric  fields  (Mozer 
et  al.  , 1977).  We  have  ignored  these  crucial  features  primarily  for  the 
sake  of  simplicity  and  for  the  saike  of  concentrating  on  the  implications  of 
pitch-angle  anisotropy,  electron  backscatter,  and  ionospheric  conditions. 
We  hope  to  remove  these  shortcomings  of  our  model  in  subsequent  works. 
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Finally,  we  take  note  of  some  very  recent  studies  by  Lennartsson 
(1976,  1977)  and  Whipple  (1977),  which  have  some  features  in  common  with 
the  present  work.  We  became  aware  of  these  studies  only  after  completion 
of  the  present  numerical  calculations,  but  we  do  not  view  the  present  work 
as  a duplication  of  their  efforts.  The  problem  of  calculating  parallel  electric 
fields  and  charged- particle  distributions  in  the  auroral  zone  is  a complicated 
one,  and  different  investigators  are  likely  to  approach  the  problem  in  some- 
what different  ways. 
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APPENDIX  A:  ACCESSIBILITY 


We  have  used  the  simple  criterion  v jjc  ? 0 in  order  to  determine 

whether  a phase- space  point  ( v^^  , v^^  ; s ) is  connected  to  a source 

point  ( '^||g*  > ^ under  the  laws  of  adiabatic  charged-particle  motion. 

Strictly  speaking,  however,  the  necessary  and  sufficient  condition  for  accessi- 
2 , ^ 

bilitv  is  that  v„  , -^0  for  all  s between  s and  s . The  equivalence  of 

' IIS 

2 

this  seemingly  more  stringent  criterion  to  the  simple  criterion  on  is 

contingent  on  the  functional  form  of  Vg(Bg),  as  will  be  shown  here. 

We  wish  to  determine  some  (preferably  minimal)  constraints  on  the 

form  of  Vg  that  will  assure  the  truth  of  the  following  statement  for  all 

2 2 2 

particles  of  interest:  if  v„  > 0 and  v.,  :jc  > 0,  then  v,,  , > 0 for 

lls  "S  lls 

* 2 
all  3 between  s and  s . The  statement  will  obviously  hold  true  if  v^^^, 

IS  a montonic  function  of  B , for  0 ^ s'  ^ c . It  will  also  hold  true  if 

s 

2 

V , is  concave  downward  when  plotted  as  a function  of  B , between  s = 0 

lls  s 

and  s'  = f . The  presence  of  q = i 1 in  (12)  would  make  it  difficult  for  both 

ions  and  electrons  to  have  a monotonic  (or  both  to  have  a concave- downward) 

2 2 2 2 

V for  all  values  of  v and  v . We  shall  therefore  try  to  make  v 

lls  lls  Is  ' lls 

monotonic  for  ions  and  concave  downward  for  electrons. 

2 

Monotonicity  of  v^^^/  for  ions  (q  = +1)  would  require  that 


(”,  ./B.) 

is  s 


- (Z/m)  )el(dV  ,/dB  ,) 
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as  can  he  shown  by  differentiating  (12)  with  respect  to  , using  a 

Cl- nte  red- dipole  model  for  the  magnetic  field  ( r , = altitude  + R_  ). 

s hj 

The  inequality  in  (Al)  is  imposed  by  the  leading  term,  which  would 

dominate  in  the  limit  of  large  . For  the  same  inequality  to  hold  in 

2 

the  limit  of  small  v,  , however,  one  must  have  dV  ,/  dB  , > 0 at 

is  s s 

least.  The  stronger  condition 

|e  (dV^,/dB^,)  > (GM^m^  /3RgBQ)(BQ/B^,)^/^  (A2) 

for  0 < s <:  i would  enable  the  electric- field  term  to  overcome  gravity, 

2 

as  is  required  for  the  unconditional  monotonicity  of  v„  , . 

' lls 

Gravity  is  negligible  for  electrons  (q  = - 1),  but  the  appearance  of 

q < 0 in  (b)  necessitates  a further  functional  constraint  on  V . We 

s 

2 

require  (see  above)  that  v for  electrons  be  concave  downward  when 

lls 

plotted  against  B^,  . In  other  words,  we  require  that 

d^(v^l^)/dB^^,  = (2/m_)  I e : (d^V^,/dBg^,)  < 0 (A3) 

for  0 S s'  I , which  is  to  say  that  V , must  be  concave 

downward  when  plotted  against  B^  . 

For  electrostatic  potentials  V^,  that  fail  to  satisfy  (A2)  and  (A3), 
there  is  a serious  danger  that  a "barrier"  in  the  "effective"  potential 
defined  (e.  g.  , Whipple,  1977)  by 
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9 . 
s 


(m  GMp,/r  ,) 


(A4) 


(m  V,  /2B  ) B , 
q 1 9 s 9 


would  prevent  the  access  of  a particle  from  s*  to  s,  despite  the  fact 
2 2 

that  both  > 0 and  > q were  satisfied  for  the  particle  in  question. 

It  IS  to  preclude  such  a danger  of  improper  mapping  that  we  require  (A2)  and 

(A3)  as  constraints  (imposed  a posteriori)  on  the  admissibility  of  "solutions" 

that  are  found  to  satisfy  (1)  for  0 S s i 1. 

We  have  carefully  checked  to  be  sure  that  our  solutions  V (B  ) 

s s 

shown  m Figures  4 and  5 satisfy  (A3),  which  is  equivalent  to  (7).  We  have 
also  verified  that  our  solutions  satisfy  (A2)  for  hydrogen  ions.  However, 
we  find  that  our  solutions  fail  to  satisfy  (A2)  at  B B for  heavier  ions 

s I 
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such  as  O . This  is  not  a serious  drawback  in  the  present  context. 

Since  we  have  considered  only  one  species  of  ionospheric  ion,  we  may 
logically  identify  the  1+  population  in  (13)  with  H and  consider  all  the 
O to  reside  beyond  s = i.  e.  , below  the  base  of  collisionless  medium. 
Alternatively,  we  may  regard  the  I ■*  population  as  consisting  of  ions  having 
a mass  equal  to  the  density- weighted  mean  of  and  m^  at  s = H . 

Either  of  the  above  alternatives  can  be  supported  by  an  appeal  to 
the  H^/O^  density  ratio  observed  at  altitudes  below  2000  km  in  high- 
latitude  trough  regions  (Hoffman  et  al.  , 1974),  but  neither  alternative  is 
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strictly  acceptable  in  the  context  of  accounting  for  the  upward  fluxes  of 
energetic  O that  are  observed  (Shelley  et  al.  , 1976)  at  high  altitudes. 

We  presume  that  a more  realistic  treatment  of  the  ionosphere  (beyond 
the  scope  of  the  present  work)  would  alleviate  this  dilemma.  Thus,  in 
a more  sophisticated  study,  one  might  treat  the  several  ionospheric 
constituents  separately  and  also  take  account  somehow  of  the  smooth 
transition  between  the  collision-dominated  lower  ionosphere  and  the 
collision- free  upper  ionosphere.  Such  a study  is  in  progress. 

The  only  O ^ ions  that  would  have  been  treated  improperly  in  the 
present  work  are  some  of  those  found  equatorward  of  the  maximum  in  the 
effective  potential  , which  is  defined  by  (A4).  Those  found  earth- 

ward of  the  maximum  in  <I>^  would  have  been  mapped  in  precise  accordance 
with  Liouville's  theorem,  except  of  course  for  the  depletion  factor  D 
given  by  (14).  It  is  not  likely  that  the  more  sophisticated  treatment  of  the 
ionosphere  described  above  will  alter  the  prsent  form  of  severely, 

so  as  to  let  O ^ ions  satisfy  (A2),  for  in  that  case  the  ionosphere  would 
rapidly  be  depleted  of  . However,  a more  realistic  treatment  of 

the  ionosphere  would  better  account  for  the  relative  distribution  of 
and  O along  the  field  line  and  might  (by  treating  the  effect  of  gravity  on 
O in  a fully  consistent  way)  enable  one  to  increase  the  convective  sweeping 
time  T that  appears  in  (14).  This  would  be  desirable. 
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APPENDIX  B:  DENSITY  MOMENT 


In  this  appendix  we  sliall  sfiow  a sample  calculation  of  the 
charge-density  moment  of  magnetospheric  electrons,  which  is  required 
for  implementing  quasi-neutrality.  With  the  use  of  the  distribution  function 
(9)  for  magnetospheric  electrons  f^^  , we  obtain  the  magnetospheric 

electron  density  moment 

n..  (s)  = / V dv  / f.  (v  ,v,  ;V  ) dv 

M-  # is  ia  I M-lis  iss  II s 

*'o  ® 

s (2.T^^_/m_)^/2  Qj^_(s)  (Bl) 

In  our  model,  the  density  of  magnetospheric  electrons  nj^  (0)  at  the  equator 
is  the  input  parameter  Nj^  ; therefore,  (Bl)  is  normalized  so  that 

~ ^'m-  ^ ^M- ^^M- ^ ’ (B2) 

where  the  function  Q,,  (0)  is  defined  as  the  limit  of  (e)  as  s -•  0 

M M ~ 

The  function  Qj^_(s)  is  given  in  terms  of  functions  related  to  error  func- 
tions and  Dawson  integrals  (Abramowitz  and  Stegun,  1974).  It  is  con- 
venient to  define 

/X 

dy  exp(-y^)  (B3) 

u 
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E^(x)  = cxp(-x^)J'  dy  exp(y^)  (B5) 

■'o 

together  with  the  simplifying  relations 

h,  = B,  - [1  - ,B6) 

>3  = (B7) 

“s  = [IB,  - BoIT_/|B„T,_)]''2 

% = [<Bi  - B^l/hj^)*'^  (B9) 

t,  - ((B,  - - B,l)  . (BlOl 

The  function  Qj^  (s)  is  given  by 


(B  /2h  )]  tr 

9 9*8 


/2)  exp[  e I 


1/2 


4 

+ e(4^)  expf.  (ypj^)2]  [ [ . 

(BID 

The  evaluations  of  the  charge  densities  are  independently 
checked  by  three  methods:  (a)  the  charge- density  integrals  for  the  mirror- 

ing and  precipitating  components  are  evaluated  separately  and  the  sum  is 
checked  against  (Bll):  (b)  if  we  set  0 equal  to  unity  in  (9),  then  there 

are  no  phase-space  restrictions  due  to  the  functions  and  the  total  charge- 

density  integral  is  easy  to  compute,  and  we  have  checked  that  the  total 
charge-density  integral  is  equal  to  the  sum  of  the  mirroring  component 
plus  twice  the  precipitating  component;  (c)  finally,  the  arguments  of  the 
0 functions  in  (9)  simplify  in  the  limit  s -•  0 ^ , and  the  evaluation  of  (Bll) 
in  the  limit  s -*  0 ^ has  been  checked  against  the  s -•  0 ^ limit  of  the 
charge-density  integral  independently  evaluated  before  the  limit  is  taken. 
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THB:  IVAN  A.  GETTING  IJ\  BORA  TOR  lES 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  (or  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  v’tal  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory;  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high- powe r gas  Lasers. 

Chemistry  and  Physics  Laboratory;  Atmospheric  reactions  and  atmos- 
pheric  optics~chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  nlumcs,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics:  quantum  electronics, 
lasers,  and  electro- optic s ; communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging,  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory;  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry,  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment,  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglowi  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields,  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 

THE  AEROSPACE  CORPORATION 
El  Segundo,  California 
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